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A bstract We continue the inv estigation of the dynamics of retrograde resonances initiated 
in M orais & Giuppond J2012h . After deriving a procedure to deduce the retrograde reso- 
nance terms from the standard expansion of the three-dimensional disturbing function, we 
specialize in the planar problem and construct surfaces of section that explore phase-space 
in the vicinity of the main retrograde resonances (2/-1, 1/-1 and 1/-2). In the case of the 1/-1 
resonance for which the standard expansion is not adequate to describe the dynamics, we 
develop a semi-analytic model based on numerical averaging of the unexpanded disturbing 
function, and show that the predicted libration modes are in agreement with the behavior 
seen in the surfaces of section. 

Keywords Resonance. Three-body problem. Surface of section. 



1 Introduction 

The discovery of extrasolar planets that orbit their host stars in the direction opposite the 
star's rot ation has renewed interest in the dynamics of retrograde motion in gravitational 
systems dTriaudet al.ll2010h . In the solar system, retrograde motion is confined to smaller 
bodies such satellites of the outer planets and long period comets. Understanding the struc- 
ture of retrograde motion and in particular re trograde resonances w ill help elucidate the 
origin and evolution of the observed systems. iGavon & Boi si J2008h performed numerical 
integrations of systems with planets moving in opposite directions and observed that a retro- 
grade resonance in two-planet systems is more stable than the equivalent prograde resonance 
confirming the idea that when two bodies orbit in different directions, encounters occur at a 
hi gher relative velocit y during a shorter time and mutual perturbations are therefore weaker. 
In lGavon et al.l 1200^) the authors obtained the 2/1 retrograde resonance disturbing function, 
and identified retrograde resonance angles. However, they were unable to find initial condi- 
tions that correspond to libration in retrograde resonance, and did not identify the theoretical 
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reason for the observed enhanced stability of retrograde resonances. InM orais & Giuppond 
d2012h (hereafter Paper I) we compared the stability of prograde and retrograde planets 
within a binary system. We observed that retrograde planets remain stable nearer to the 
secondary star than prograde planets. We showed that instability is caused by single mean 
motion resonances (MMRs) and the possible overlap of adjacent pairs. We used a standard 
expansion of the disturbing function for the planar circular restricted three-body problem 
(CR3BP) to obtain the retrograde resonance terms, and we explained how these terms show 
why retrograde resonances are more stable than prograde resonances as the magnitude of 
the p/q resonance terms is proportional to a power of the eccentricity, which at the lowest 
order, is e p+q in retrograde case and e p ~ q in prograde case. 

In this paper, we continue our investigation of retrograde motion in the three-body prob- 
lem by systematically studying the structure of the phase space near the main retrograde 
resonances. We concentrate on the planar problem and examine in detail motion near the 
2/1, 1/1 and 1/2 retrograde resonances. In section 2, we explain how to obtain the retro- 
grade resonance terms from an expansion in Laplace coefficients of the three-dimensional 
disturbing function. We also show that, in the planar problem, these coincide with the ret- 
rograde resonance terms obtained in Paper I. Since Laplace coefficients diverge when the 
semi major axes' ratio is close to unity, in section 3, we develop a semi-analytical model 
for the co-orbital 1/1 retrograde resonance based on numerical averaging of the unexpanded 
disturbing function. In section 4, we present our numerical approach and results and de- 
scribe how retrograde motion phase-space is structured and where stable motion is possible. 
Section 5 contains a discussion of our results. 



2 Differences between prograde and retrograde resonance 

The encounter of two bodies in a retrograde configuration (orbiting in opposite directions) 
occurs at a higher relative velocity during a shorter time than in a prograde configuration. 
This implies that mutual perturbations are weaker for retrograde MMRs. In Paper I, we 
studied retrograde MMRs analytically in the context of the planar circular restricted three- 
body problem and compared their relative strength and stability to prograde resonances. 
Here, we will explain how we obtain the slow terms of the disturbing function for a p/q 
retrograde MMR in the three-dimensional CR3BP. 



2.1 Disturbing function 

Consider a test particle that moves under the gravitational effect of a binary composed of a 
primary with mass Af* and a secondary with mass m <C M* . The motion of m with respect 
to M* is a circular orbit of radius a' = 1 and longitude angle A'. The reference plane 
is defined by the binary's orbit. The test particle's osculating Keplerian orbit with respect 
to M* has semi-major axis o, eccentricity e, inclination /, true anomaly /, argument of 
pericentre u, and longitude of ascending node Q. The disturbing function reads: 

R = Gm{ l/A - r cosip) (1) 

where r is the radius of test particle, ip is the angle between the radius vectors of binary and 
test particle, A 2 = 1 + r 2 — 2 r cos tp, and 



cosf/' = cos(f2 — A') cos(oj + /) — sm{fl — A') sin(aj + f) cos I. (2) 



Retrograde resonance in the planar three-body problem 



3 



The first term of R (direct perturbation) is the gravitational force from the mass m on the 
test particle whereas the second term (indirect perturbation) comes from the reflex motion 
of the star under the influence of the mass m as the standard coordinate system is chosen to 
be centered on the star. 

The classic series of the disturbing function is expanded in powers of sin 2 (//2). This 
is adequate for nearly coplanar prograde motion since sin 2 (1/2) w but not for nearly 
coplanar retrograde motion since sin 2 (//2) w 1. Therefore, for nearly coplanar orbits, we 
define /3 < 1 such that I = f3 or / = 180° — /3, respectively, for prograde or retrograde 
motion. We may therefore write: 

cosi/> = (1 - s 2 ) cos(/ + oj ± (f2 - A')) + s 2 cos(/ + uj =p {O - A')), (3) 

where s 2 = sin 2 (/3/2) <C 1, and the ± sign applies to prograde or retrograde cases, respec- 
tively. 

Next, we write A 2 = 1 + r 2 - 2 r cos(/ + u ± (Q - A')) - 2r&, where if - is defined 

as: 

tf = cos-0 - cos(/ + uj± (Q - A')) = 2s 2 sin(±(^ - A'))sin(w + /). (4) 

Expanding the direct perturbation term A -1 in the vicinity of = (as s 2 <C 1), we may 
write: 

A z. (l!)2 ^ 2 ; ai^ Pj 



where A 2 = 1 + r 2 - 2 r cos(/ + u ± {Q - \') 
r j a — 1 = 0(e) and exp 

1 / , . \ ^ 1 k k d k 



Finally, defining e = r/a — 1 = O(e) and expanding A ^ 2l+1 ' around e = 0: 



l + J2h kak 7L (6) 



A 2 t+1 \ it 1 ! kl dka J p2 

with a = a/a' , and 

_L_ = (l + Q 2 -2acos(/ + ^±(f2-A')))- (l+1/2) , 

= E ^+l/ 2 («) cos (j(/ + w ± (fi - A'))), (7) 
j 

where fr? | /? (q) are Laplace coe fficients. For a < 1, they may be expanded as convergent 
series in a fellis & Murravll2000l) . 

Examination of the angles in the expression of *P (O and that of p Q show that there are 
two ways with which the expansion of the disturbing function for retrograde motion can be 
obtained from the expansion of the disturbing function for prograde motion. The two ways 
are equivalent as they depend on whether one chooses to invert the motion of the inner body 
or that of the outer on43- Inverting the motion of the outer body gives the first transformation: 
I* = 180° - /, A'* = -A', uj* = uj - 7r and J?* = -J? - 7r. In this case, the longitude 
of pericentre vj = uj + fl is transformed into vj* = uj — fi and the mean longitude A = 



1 A retrograde orbit with inclination / > 90° can be obtained from a prograde orbit with inclination 
180° — / by inverting the direction of motion which implies a swap between ascending and descending 
nodes. 
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M + lo + f2 into A* = M + u — Q where M is the mean anomaly. Equivalently, this means 
applying the generating function F\ = -X'A'* + (A + 2z)A* + (g - 2z)T* - (z - tt)Z* 
to the usual Poincare action-angle variables. Inverting the motion of the inner body gives the 
second possible transformation: I* = 180° — I, f* = — /, uj* = n — uj and fl* = ir + fl 
that may be obtained with the generating function F 2 = X'A'* - (\+2z)A* - (g - 2z)T* + 
(z — ir)Z* . We note however that these two transformations are passive in that they allow 
us only to obtain the expression of the resonant arguments. Once the arguments are obtained 
formally, the assumption that A > and A' > always holds. This is in contrast to the 
approach adopted in Paper I where an active transformation was used to study the planar 
dynamical problem by choosing explicitly from the outset the convention A > and A' < 0. 

2.2 Resonant terms 

Now that we have shown how the expansion of the disturbing function for J = 180° — j3 
is obt ained from that with i" = /3, we may use the literal expansion of lEllis &Murravl 
d2000h valid for prograde motion and transform the relevant resonance terms to describe the 
corresponding retrograde resonance. 

The 2/1 retrograde resonance terms are of type e 3 cos( A — 2 A' — 3 zu) [term 4D3. 1 with 
j = 2], and e s 2 cos(A - 2 A' - vj + 2 fl) [term 4D3.5 with j = 2]. The 1/2 retrograde 
resonance has direct and indirect terms of type e 3 cos(2 A — A' — 3 vj) [terms 4D3.4 with j = 

2 and 413.6], and e s 2 cos(2 A - A' - zu + 2 Q) [terms 4D3.10 with j = 2 and 413.13]. The 
1/1 retrograde resonance has direct and indirect terms. These are of type e 2 cos(A — A' — 2 vj) 
[terms 4D2.1 with j = 1 and 4E2.2 ], and s 2 cos(A - A' + 2 Q) [terms 4D2.4 with j = 1 
and 4E2.6]. However, the 1/1 resonance direct terms cannot be obtained from the literal 
expansion of the disturbing function (since Laplace coefficients diverge when a — > 1). We 
develop a semi-analytic model for the co-orbital resonance in the next section. 

A similar analysis for any p/q retrograde resonance shows that there are resonant termsQ 

e P+q-2 k g 2 k cos ( 9A _ pX > _ ( p + q _ 2 + 2ktt). (8) 

with k = 0,1, 2, ... andp + q > 2k. 

li /3 <ti I then s <C 1, hence the term with k = is dominant. Since we restrict our 
study to planar retrograde resonance (s = 0), only the term e p+? cos <j> remains, where 

4> = <?A - pX' - (p + q)zu. (9) 

We thus recover the retrograde resonant angle from Paper I (with the expected change of sign 
for the term in A'). Following Paper I, we use the notation p/-q resonance when referring to 
a p/q retrograde resonance. 

3 A model for co-orbital resonance 

Consider a test particle in the co-orbital region of the secondary (| a— 1 1 <C 1). The disturbing 
function may be expressed using the natural angles: the fast epicyclic motion represented by 
the mean longitude A of the particle and the guiding centre phase represented by the relative 
mean longitude r = X — X' . To obtain the resonance hamiltonian, the disturbing function, 



2 D'Alembert rule is not obeyed because the assumption A' > and A > with vj = uj — Q implies 
that we measure angles for the test particle in the opposite direction of the binary's motion. 
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R (Eq. 1), is averaged with respect to the fast angle A. The corresponding function is the 
ponderomotive potential S = (R) used in our previous work on the co-orbital resonance 
dNamouml 1 19991 ; iNamouni et al.L Il999h . When the relative longitude r is introduced, we 
may write: 

cosip = - (1 + cos/)cos(/ -M + t) + - (1 - cos/)cos(/ + M -t + 2lj), (10) 
and the ponderomotive potential is given as: 

where r = a(l — e 2 )/(l + ecos/), and the average over the mean anomaly M has been 
replaced by an average over / using the conservation of angular momentum. The mean 
anomaly M is related to the true longitude by the eccentric anomaly tan(_E/2) = (1 — 
e) 1/2 (l + e)" 1/2 tan(//2) and Kepler's equation M = E — esinE. 

The expansion-free expression of ip gives the natural resonant angles for planar motion. 
For prograde motion (cos I = 1), libration occurs around d> = r whereas for retrograde 
motion (cos / = —1), libration occurs around (f> = r — 2ojJ. FigureQ](top row) shows the 
shape of the potential S as a function of the resonant argument </> = r — 2uj for planar 
retrograde motion where a — 1 = 0. 01. At low eccentricity, libration occurs only around 
4> = 180° and the potential is quite shallow. This explains why in the next section we ob- 
serve that low eccentricity libration orbits for relatively large mass ratios (e.g. /i = 0.01) 
are difficult to set up as the larger the mass ratio the stronger the mutual perturbations, the 
more destructive the close encounters. We shall show that such librations are quite stable 
at smaller mass ratios. As the eccentricity is increased, the collision boundary appears and 
librations may occur around or 180°. The extent of the libration amplitude depend on 
eccentricity, for 0.1 < e < 0.7, libration around zero has the largest amplitude. Libration 
around 180° regain some importance as e approaches unity. As the ponderomotive potential 
S is derived for three-dimensional orbits, it is instructive to see how a the introduction of a 
small inclination modifies the dynamics as realistic orbits in the planetary three-body prob- 
lem never lie exactly on the same plane. Moreover, in the co-orbital resonance, inclination 
is know to mitigate collision al encounters and facilitate stable orbital transitions dNamounil 
ll999l : lNamouni et al.Lll999l) . Figure [T] (middle and bottom rows) shows how S is modified 
when the retrograde orbits have a mutual inclination of 10° and ui is set to zero. As ex- 
pected, collision singularities are absent. There appears a bifurcation near e = 0.161 where 
librations around and 180° have comparable amplitudes and may be associated with the 
same energy level. The remaining features of the planar problem are present: the potential's 
shallowness for low eccentricity and the dominance of libration around zero for more eccen- 
tric orbits. We remark that the potential's amplitude and bifurcation modes depend on the 
relative semi-major axis. Figure [T]only illustrates the similarities and differences with the 
planar problem. We also note that in the fully three dimensional problem, the time evolution 
of the argument of pericenter modifies the potential's shape and equilibria whereas for pla- 
nar orbits, the potential depends on the combined phase <f> = r — 2uj. We shall present the 
study of the three-dimensional retrograde co-orbital resonance elsewhere. 



3 Here, we define A = M + zu with zu = oj + fl. If we define A = M + zu with zu = ui — Q the 
retrograde resonant angle is tj> = A — A' — 2 zu in agreement with the conclusions of the previous section. 
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4 Retrograde resonance phase space 

Poincare surfaces of section are useful tools to study the phase space structure in the three- 
body problem. In what follows, we define how we set up surfaces of section for the phase 
space of the 2/-1, 1/-1 and 1/-2 resonances. We then examine the types of orbits involved, as 
well as they potential stability. 



4.1 Surface of section construction 
In the barycentric rotating frame, the planar circular restricted three-bod y problem has two 



degre es of freedom (x, y) and one integral of motion, the Jacobi constant dMurrav & Dermotn, 
Il999ri : 

n 2 i 2 , -2n , 2(l-/i) 2/1 

C = x + y - (x +y)-\ 1 , (12) 

where n < 1 is the mass ratio of the secondary and the primary, r\ = {x — /i) + y 2 and 
T2 = {x— l + fi) 2 + y 2 . Orbits therefore lie on a 3D subspace C(x,y,x,y) = C embedded 
in the 4D phase space. Points of an orbit that intersects a given surface, e.g. y = 0, in a given 
direction, e.g. y > 0, lie on a 2D surface of sec tion (x, x). An order k res onance corresponds 
to a set of k islands on the surface of section dWinter & Murravlll997allbh . Here, we prefer 
to define the surface of section by x = 0, allowing us to follow orbit intersections in the 
(x,y)-plane. 

We choose a mass ratio fi = 0.01 that is small enough for perturbation theory to apply 
and keplerian osculating elements to be used. These elements (o, e, w) vary on a longer 
scale than the orbital period. Setting xo = 0, we have: 



^ + yo 2 + rrt+ h f\ M -c. (is) 

xo + Ml \l-(x +fj,)\ 

The transformation between the barycentric frame and astro centric frame (centered on the 
primary) at t = is given as: 

xi = xq + fj, xi = xo — yi (14) 

yi = yo yi = yo + xi (15) 

When ±o=0 and yo = we have x\ = and y\ = yo + x\. Hence, the necessary 
and sufficient conditions for setting up prograde orbits are xi > and yo > or x\ < 
and yo < (since y\ > if x\ > and yi < if xi < ). However, the necessary 
and sufficient conditions for retrograde orbits are xi > and yo < —xi or xi < and 
yo > —Xi (since > yo > —xi and x\ > implies yi > 0, whereas < yo < —x\ 
and xi < implies y\ > 0, i.e. these orbits are prograde). By substituting the critical value 
yo = —x\ in Eq. d X31) . we obtain a limit on the Jacobi constant 

C< 2 , (1 ~^ + | ^— -M(M + 2x ), (16) 

\xq + fi\ \xo — 1 + n\ 

such that, within this range of C, Xi > and yo < or xi < and yo > are necessary 
and sufficient conditions for retrograde orbits. When \xi \ < 3, this limit is C < 0.7. 

We construct surfaces of section defined by x = and y x yo > so that the initial 
condition lies on the surface of section. We vary xo between -3 and 3 with increments of 
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0. 05, and < yo < 0.8 with increment of 0.1. The Jacobi constant —1.9 < C < 0.7 is 
incremented by 0.1. We use Eq. d!3fc to obtain yo and choose yo > if x\ < and yo < 
if X\ > 0. This ensures that the initial conditions always correspond to retrograde orbits (as 
C < 0.7). The integration time is 5 x 10 binary periods. These surfaces of section can be 
seen as electronic supplementary material (a selection will be discussed in Sect. 3.3). 

We show below the level curves of constant C in (a, e) space for initial conditions at 
conjunction i.e. x\ > and y = (Fig.[2](a) and (b)) and for initial conditions at opposition 

1. e. xi < and y = (Fig.[2](c) and (d)) . The chosen range of C spans semi-major axes 
between 0.5 and 1.5 to include the 2/-1, 1/-1 and 1/-2 resonance regions. 



4.2 Resonant angles 

Test particles are started with A = A' = (conjunction: y = and x + /i > 0) or A = 180° 
and A' = (opposition: y = and x + fx < 0). The points on the surface of section 
(i = 0) with y = correspond to the osculating orbits' pericenter (A = zu) or apocenter 
(A = vj + 180°), depending on the x and C values. 

Starting in conjunction and at pericenter (apocenter) corresponds to the resonant angle 
cj) = q A-pA' - {p + q)vj = {[p + q] x 180°). In the latter case, = (180°) if p + q is 
even (odd). Starting at opposition and at pericenter (apocenter), <f> = —p x 180° (q x 180°). 
Hence for opposition at pericenter, cf> = (180°) if p is even (odd) and at apocenter = 
(180°) if q is even (odd). Therefore, even order resonant angles (such as that of the 1/-1 
resonance) librate around for initial conditions at conjunction, and librate around 180° for 
initial conditions at opposition. Odd order resonant angles with p even (such as that of the 
2/-1 resonance) librate around 180° for initial conditions at apocenter, and librate around 
for initial conditions at pericentre. Odd order resonant angles with an odd p (such as that of 
the 1/-2 resonance) librate around 180° for initial conditions at conjunction and apocenter), 
or opposition and pericenter, and librate around for initial conditions at conjunction and 
pericentre, or opposition and apocenter. 



4.3 Results of numerical integrations 
4.3.1 Resonant configurations 

We examine a selection of orbits in the 2/-1, 1/-1 and 1/-2 resonances shown in the frame 
rotating with the binary. 

Fig. <[3]> shows 2/-1 resonant orbits. The top left panel shows an orbit with C = +0.6 
that starts at conjunction or opposition (y = 0), and pericenter (mode A). The top right 
panel shows an orbit with C = +0.6 that starts at conjunction or opposition (y = 0) and 
apocenter (mode B). Resonant libration occurs around in mode A and around 180° in 
mode B (Fig. [7). The low left panel shows an orbit with C = +0.3 that starts at conjunction 
or opposition (y = 0) and pericenter (mode A). The low right panel shows an orbit with 
C = +0.3 that starts at conjunction or opposition (y = 0) and apocenter (mode B). The 
latter orbit is very close to collision with the secondary. 

Fig. © shows 1/-1 resonant orbits. The top left panel shows an orbit with C = +0.6 
that starts at x\ > and y = or xi <0 and y ^ (mode I). Resonant libration occurs 
around and disruptive close encounters are avoided despite the high eccentricity (Fig. [8] 
left: mode I). The top right panel shows an orbit with C = +0.6 that starts at x\ < and 
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y = (mode II). The 1/-1 resonant angle librates around 180° (Fig.[8]right: mode II) but the 
orbit is close to collision and becomes unstable when C < 0.6. Both orbits are described by 
the equilibria of the ponderomotive potential S in Fig. (1, upper row, rightmost panel). 

The mid left panel of Fig. © shows an orbit with C = —0.9 and moderate eccentricity 
that starts at x\ > and y = or x\ < and y ^ (mode I). The resonant angle 
librates <j> around 0. The mid right panel of Fig. © shows an orbit with C = —1.1 that is 
nearly circular and has initially x\ < and y = (mode III). This crossing orbit is very 
close to collision and we expect resonant libration around 180° (Fig. [9] right: mode III). 
What happens is an interesting behavior best seen if we integrate a similar orbit for smaller 
mass ratios thus reducing the jitter due to close encounters. In Fig. ([5] left panel), we plot 
in the ((/>, e)-plane, a similar orbit but with [i = 10 -4 . The resonant argument alternates 
periodically between libration and circulation in a state that is stable over long time scales. 
Observing libration around (f> = 180° requires a finer search which becomes easier as the 
mass ratio is decreased (Fig. [5] right panel). The low left panel of Fig. © shows an orbit 
with C = —1.2 that has small eccentricity and starts at x\ > and y = or x\ < and 
y (mode I). This crossing orbit is also close to collision and the resonant angle cf> librates 
around (Fig.|9]left: mode I). The low right panel of Fig. (|4j shows an orbit with C = —1.2 
that is nearly circular and starts at a;i > Oorxi < and y = for C = —1.2. This is a 
non-crossing orbit just exterior to the secondary's orbit and the resonant angle circulates. 

Fig. © shows 1/-2 resonant orbits. The top left panel shows an orbit with C = —1.5 
that starts at conjunction and pericenter, or opposition and apocenter (mode A). The top 
right panel shows an orbit with C = —1.5 that starts at conjunction and apocenter, or 
opposition and pericenter (mode B). Resonant libration occurs around in mode A and 
around 180° in mode B (Fig. 1 10b. When C = —1.5 both mode A and mode B orbits are 
close to collision with the secondary. The low left panel shows an orbit with C = —1.2 that 
starts at conjunction and pericenter, or opposition and apocenter (mode A). This is the only 
stable configuration when C > —1.3. The low right panel shows an orbit with C = —1.8 
that starts at conjunction and apocenter, or opposition and pericenter (mode B). This is the 
only stable configuration when C = —1.8. 

4.3.2 Surfaces of section 

As seen in Figs.[3]S|6] an order k resonant orbit intersects the section x = at 2k different 
points. However, owing to the constraint on the sign of y, we only see k of these intersections 
on the surface of section {x\, y). Therefore, in the surfaces of section (Fig. [TT1 a set of k 
intersections on the left hand side (x± < 0) usually represents the same configuration (in 
the synodic frame) as a set of k intersections on the right hand side (x± > 0). 

Fig. II II shows a selection of surfaces of section with C = 0.6, C = 0.3, C = 0.0, 
C = -0.9, C = -1.1, C = -1.2, C = -1.5 and C = -1.8. The full set of surfaces 
of section (C between 0.7 and —1.9 at steps 0.1) can be seen as electronic supplementary 
material (ESM). At these values of C many initial conditions correspond to crossing orbits 
(see Fig. 2) which can only be stable in resonance. Non-crossing small eccentricity orbits 
exist in the regions marked in green (left and right on Fig. [PH . 

When C = 0.6, we see nearly circular non-crossing orbits in the vicinity of the 3/-1 
resonance (green) for initial conditions at opposition (left) or conjunction (right). Collision 
with the secondary occurs between the 3/-1 and 2/-1 resonances. We see islands of libration 
in 2/-1 resonance. The 2/-1 resonant orbits that start at pericenter (magenta on left and right) 
correspond to the same configuration in the synodic frame (Fig.[3]top left panel: mode A) 
where the resonant angle librates around (Fig. |7]left). The 2/-1 resonant orbits that start at 
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apocenter (blue on left and right) also correspond to the same configuration in the synodic 
frame (Fig.|3]top right panel: mode B) where the resonant angle librates around 180° (Fig. [7] 
right). 

When C = 0.6, we also see islands of libration in 1/-1 resonance. The 1/-1 resonant 
orbits that start at conjunction (black on right) or with x\ < and y ^ (black on left) 
correspond to the same configuration in the synodic frame (Fig.|4]top left: mode I) where the 
resonant angle librates around (Fig.|8]left). The 1/-1 resonant orbit that starts at opposition 
(red) is close to collision (Fig. |4]top right: mode II) and the resonant angle librates around 
180° (Fig.[8]right). The 1/-1 resonance occurs at very high values of the eccentricity and the 
resonance center has a = 1. 

When C = 0.3, we see orbits the vicinity of the 2/-1 resonance for initial conditions 
at opposition (left) or conjunction (right). The 2/-1 resonant angle can circulate for non- 
crossing orbits (green on left and right), it can librate in mode A (Fig. [3] low left) for initial 
conditions at pericenter (magenta on left and right), or it can librate in mode B (Fig. [3] low 
right) for initial conditions at apocenter (blue on left and right). Collision with the secondary 
occurs just outside the 2/-1 resonance separatrix. The 1/-1 resonant angle librates around 
(mode I) for initial condition at conjunction (black on right) or with x\ < and y ^ 
(black on left). There are also islands of libration in the 1/-2 resonance for initial condition 
at conjunction / pericenter (magenta on right) or opposition / apocenter (magenta on left). 
These correspond to the same configuration in the synodic frame where the resonant angle 
librates around (Fig.|6]mode A). 

When C = 0.0, we see nearly circular orbits in the vicinity of the 2/-1 resonance (green) 
for initial conditions at opposition (left) or conjunction (right). The 2/-1 resonant angle can 
only circulate and all these orbits are non-crossing. Collision with the secondary occurs 
between the 2/-1 and 3/-1 resonances. The 1/-1 resonant angle librates around (mode I) 
for initial condition at conjunction (black on right) or with x± < and y ^ (black on 
left). The 1/-2 resonant angle librates around (mode A) for initial condition at conjunction 
/ pericenter (magenta on right) or opposition /apocenter (magenta on left). 

When C = —0.9, all initial conditions correspond to crossing orbits hence they are only 
stable in resonance. The 1/-1 resonant angle librates around (mode I) for initial condition 
at conjunction (black on right) or with x\ < and y ^ (black on left). Collision with 
the secondary occurs in the 1/-1 resonance region. The 1/-2 resonant angle librates around 
(mode A) for initial condition at conjunction / pericenter (magenta on right) or opposition 
/apocenter (magenta on left). 

When C = — 1.1, nearly circular orbits in the 1/-1 resonance starting at opposition (red 
on left) correspond to libration around 180° (Fig. [4] mid right: mode III) although that is 
not very clear from the behavior of the resonant angle (Fig. |9] right: mode III) due to the 
effect of repeated very close encounters. These are crossing orbits which are very close 
to the collision boundary. The nearly circular orbit starting at conjunction (red on right) 
correspond to the 1/-1 resonance separatrix. The 1/-1 resonant angle can also librate around 
(mode I) for initial conditions at conjunction (black on left) or with x± < and y ^ 
(black on right). These are also crossing orbits close to the collision boundary. The 1/-2 
resonant angle librates around (mode A) for initial condition at conjunction / pericenter 
(magenta on right) or opposition /apocenter (magenta on left). 

When C = —1.2, we see nearly circular orbits in the vicinity of the 1/-1 resonance 
(green) for initial conditions at opposition (left) or conjunction (right). The 1/-1 resonant 
angle circulates and the orbits are just exterior to the secondary's orbit thus very close to 
the collision boundary (Fig.|4]low right: mode III). The 1/-1 resonant angle can also librate 
around (Fig. [9]left and Fig.|4]low left: mode I) for initial conditions at conjunction (black 
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on right) or with xi < and y ^ (black on left). These are crossing orbits close to the 
collision boundary. The 1/-1 resonance is no longer possible when C = —1.3 (see ESM). 
When C = —1.2, the 1/-2 resonant angle librates around (mode A) for initial conditions 
at conjunction / pericenter (magenta on right) or opposition /apocenter (magenta on left). 
When C = —1.3 (see ESM) the 1/-2 resonant angle can also librate around 180° (mode B) 
for initial conditions at apocenter / conjunction and pericenter / opposition. 

When C = —1.5, we see nearly circular orbits in the vicinity of the 2/-3 resonance 
(green) for initial conditions at opposition (left) or conjunction (right). Collision with the 
secondary occurs in the 21-3 resonance region. There are islands of libration in 1/-2 reso- 
nance. The 1/-2 resonant orbits that start at pericenter / conjunction (magenta on right) and 
apocenter / opposition (magenta on left) correspond to the same configuration in the synodic 
frame (Fig.|6]top left panel: mode A) where resonant angle librates around (Fig. [lO]left) . 
The 1/-2 resonant orbits that start at apocenter / conjunction (blue on right) and pericenter / 
opposition (blue on left) correspond to the same configuration in the synodic frame (Fig. [6] 
top right panel: mode B) where the resonant angle librates around 180° (Fig. [TOl right). 

When C = —1.8, there are nearly circular orbits in the vicinity of the 1/-2 resonance 
(green) for initial conditions at opposition (left) or conjunction (right). The 1/-2 resonant 
angle can circulate or it can librate around 180° (mode B) for initial condition at opposi- 
tion / pericenter (blue on left) or conjunction / apocenter (blue on right). Collision with the 
secondary occurs just outside the 1/-2 resonance separatrix. The 1/-2 resonant angle can no 
longer librate around (mode A). 



4.4 Analytic model for 2/-1 and 1/-2 resonances 

The structure of the 2/-1 and 1/-2 resonances, at low to moderate ec centricities, can be 
descri bed by the analytic model for 3rd order resonance presented in iMurrav & Dermotj 
d 19991) and described in Paper I. The Hamiltonian is: 

H = ^(X 2 + Y 2 ) + i(X 2 +Y 2 ) 2 t2X(X 2 - 3Y 2 ) (17) 

where 5 measures the proximity to exact resonance, =p applies to 2/-1 or 1/-2 resonances, 
X = Rcos((f>/3), Y = i?sin((ft/ 3) and R is a scaling facto r specific for each resonance 



and dependent on the mass ratio /i dMurrav & Dermoid 1 19991) . For the 21- 1 resonance 

2 T 

A — 2 A — 3 vj and R = 2 3 7 e, where 7 = (3/i fs2) and fs2 = 0.402 is the amplitude 
of term 4D3.1 (j = 2) when a = 0.623 . For the 1/-2 resonance cf> = 2 \ — A' — 3ro 
and R = 2 2 7 e where 7 = (3fi (/ 85 - 0.5/q 2 ))" 1 where / 85 - 0.5/a 2 = 0.533 is the 
combined amplitude of the terms 4D3.4 (j = 2) and 413.6 when a = 0.623. 

A resonant orbit corresponds to a set of 3 stable equilibrium points of the Hamiltonian 
(Eq.fTTll. In Paper I we showed curves of constant Hamiltonian for several values of the pa- 
rameter 5. When 5 = (exact resonance) there is a bifurcation at the origin, and the 3 stable 
equilibrium points have R = 6 and (f> = 0, ±27r/3 for 2/-1 resonance and <j> = n, ±7r/3 for 
1/-2 resonance. Applying the scaling we see that, when fi = 0.01, the equilibrium points 
at the 2/-1 resonance have e = 0.05, while the equilibrium points at the 1/-2 resonance 
have e = 0.03. In Fig. [12] we show real trajectories in the vicinity of the 2/-1 resonance 
and in 1/-2 resonance obtained by numerical integration of the equations of motion with 
fi = 0.01 at Jacobi constant values (C) close to bifurcation at the origin. There is very good 
agreement with the analytic model for the 2/-1 resonance (Fig |12ab the scaling is correct 
and the topology in (X, Y) coordinates is similar to the topology of the surface of section 
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(same equilibrium points and separatrices ). However, for the 1/-2 resonance the scaling is 
still approximately correct but a set of 3 stable equilibrium points in the surface of section 
correspond to the large amplitude libration orbit in Fig J12bl possibly due to the vicinity of 
the separatrix at e = 0. 



5 Discussion 

This article is the continuation of our work on retrograde resonances (Paper I). We identified 
the transformation that must be applied to the standard expansion of the three-dimensional 
disturbing function in order to obtain the relevant resonance terms for retrograde motion and 
analyze them quantitatively. 

We explored the phase-space near the retrograde resonances 1/1, 2/1, and 1/2, by con- 
structing surfaces of section for the planar three-body problem using a mass ratio 0.01. The 
1/-1 resonant term amplitude is of order e 2 and the 2/-1 and 1/-2 amplitudes are of order e 3 , 
hence these are the strongest retrograde resonances. 

We saw that for low eccentricity non-crossing orbits, libration in the 2/-1 resonance is 
around (corresponds to starting at conjunction or opposition and pericenter) whereas li- 
bration in the 1/-2 resonance is around 180° (corresponds to starting at conjunction and 
apocenter, or opposition and pericenter). These are the most stable configurations for non- 
crossing resonant orbits, since they ensure that closest approach with the secondary occurs 
always at pericenter for the 2/-1 resonance, and always at apocenter for the 1/-2 resonance. 
The behavior of low eccentricity 2/-1 and 1/-2 resonant orbits is in reasonable agreement 
with an analytic model for 3rd order resonances based on the literal expansion of the dis- 
turbing function. However, this analytic model is not valid for resonant crossing orbits. We 
observed that moderate to large eccentricity crossing orbits in the 2/-1 resonance or mod- 
erate eccentricity crossing orbits in the 1/-2 resonance can librate around or 180°, while 
large eccentricity crossing orbits in the 1/-2 resonance can librate only around 0. 

Recalling that the literal expansion of the disturbing function is not adequate for co- 
orbital motion, we developed a semi-analytic model for the 1/1 retrograde resonance based 
on numerical averaging of the full disturbing function and valid for large eccentricity and 
inclination. We saw that this model correctly explains the 1/-1 resonant modes, namely li- 
bration around and around 180°. Whereas 1/-1 resonant libration around is quite stable 
and occurs between e w 0.2 and e w 1, libration around 180° occurs only at e « and 
e « 1, hence it is located close to the collision separatrix with the secondary. The disruptive 
effect of these close encounters implies that 1/-1 libration around 180° is easier to setup for 
smaller mass ratios (< 10 -4 ). We also expect that transitions between different 1/1 retro- 
grade resonant modes are possible in the thr ee-dimensional problem, in analogy wi th what 
is described for prograde 1/1 resonant orbits (lNamouniUl999l : lNamouni et all 1999h . 
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Fig. 1 Ponderomotive potential S of the 1/-1 resonance as a function of the resonant angle. The relative 
semi-major axis is 0.01. The upper row deal illustrates the planar potential whereas the middle and bottom 
rows, three dimensional orbits with small inclination. For the latter, r — 2u) is no longer the only possible 
resonant argument and S is plotted as a function of r for a fixed to. 
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Fig. 2 Level curves of C at values (from left to right) 0.6,0.3,0,-0.3,-0.6,-0.9,-1.2,-1.5,-1.8. Initial condition 
at conjunction and pericenter (a), conjunction and apocenter (b), opposition and pericenter (c), opposition 
and apocenter (d). The magenta and blue lines locate collision with secondary at pericenter or apocenter. The 
red vertical lines show location of 2/-1, 1/-1 and 1/-2 resonances. 
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Fig. 4 Orbits in 1/-1 resonance seen in synodic frame: C = 0.6 mode I (top left) and mode II (top right); 
C = -0.9 mode I (mid left); C = —1.1 mode III libration (mid right) C = -1.2 mode I (low left) and 
mode III circulation (low right). A unit radius circle in cyan helps identify the crossing orbits and non-crossing 
orbits. 
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Fig. 6 Orbits in 1/-2 resonance seen in synodic frame: C = —1.5 mode A (top left) and mode B (top right); 
C = -1.2 mode A (low left) and C = -1.8 mode B (low right). 
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Fig. 7 Orbits in 2/-1 resonance: C = 0.6 mode A (left) and mode B (right). 
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Fig. 8 Orbits in 1/-1 resonance: C = 0.6 mode I (left) and mode II (right). 



18 



M. H. M. Morais, F. Namouni 



a=1. 1513246716779. e=0.19814742643455, 1=180°, M=89.68°, co=72.95° a=1 .0692580589982, e=1.93984425248879E-002, 1-180°, M=180°, to=0 



200 400 600 

time 




1000 200 400 600 800 1000 



200 400 600 800 1000 200 400 600 800 1000 



360 
270 
E 180 



360 
270 
180 



■tri 



200 400 600 800 1000 200 400 600 800 1000 

time time 



360 
270 



360 
270 



200 400 600 

time 



1000 200 400 600 800 1000 



Fig. 9 Orbits in 1/-1 resonance: C = —1.2 mode I (left) and C = -1.1 mode III (right). 
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Fig. 10 Orbits in 1/-2 resonance: C = —1.5 mode A (left) and mode B (right). 
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Fig. 11 Surfaces of section for selected values of C. Primary is at (0, 0) and secondary is at (1, 0). 
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Fig. 11 Surfaces of section for selected values of C. Primary is at (0, 0) and secondary is at (1, 0). 
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Fig. 12 (a) trajectories in the vicinity of 2/-1 resonance obtained by numerical integration of the equations 
of motion with /i = 0.01 at C = 0.051: near exact resonance (black), separatrix (red and green) and outer 
circulation (magenta), (b) trajectories in the vicinity of 1/-2 resonance obtained by numerical integration of 
the equations of motion with fi = 0.01 at C = —1.886: near exact resonance in surface of section (red). 



